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ABSTRACT 

Recent Chandra observations have revealed a large population of faint X-ray point sources in the 
Galactic Centre. The observed population consists of ;> 2000 faint sources in the luminosity range ~ 
10 31 -10 33 erg s _1 . The majority of these sources (70%) are described by hard spectra, while the rest are 
rather soft. The nature of these sources still remains unknown. Belczynski & Taam (2004) demonstrated 
that X-ray binaries with neutron star or black hole accretors may account for most of the soft sources, 
but are not numerous enough to account for the observed number and X-ray properties of the faint hard 
sources. A population synthesis calculation of the Galactic Centre region has been carried out. Our 
results indicate that the numbers and X-ray luminosities of intermediate polars are consistent with the 
observed faint hard Galactic Centre population. 

Subject headings: Galaxy: center — X-rays: binaries — stars: white dwarfs 



1. INTRODUCTION 

A Chandra X-ray survey of the Galactic Centre (GC, 
Wang, Gotthelf & Lang 2002) first revealed the pres- 
ence of ~ 1000 spectrally hard X-ray sources (2-10 keV) 
with luminosities L x <J 10 35 erg s _1 . Pfahl et al. 
(2002) have claimed that wind-fed neutron stars (NS) with 
intermediate- and high-mass companions are responsible 
for a significant fraction of the hard sources in the Wang et 
al. (2002) survey. A deeper Chandra survey of the nuclear 
bulge region (Muno et al. 2003; hereafter MM03) revealed 
over 2000 X-ray point sources with X-ray luminosities ~ 
10 31 -10 33 erg s" 1 . The majority of these sources (1427) 
are described by hard spectra (with a photon index of an 
absorbed power law F < 1), while the rest (652) are char- 
acterized by softer spectra (r > l). 1 Belczynski & Taam 
(2004; hereafter BT04) have studied the entire population 
of X-ray binaries with NS and black hole (BH) accretors in 
the context of the MM03 survey. It was demonstrated that 
neither wind- fed (low-, intermediate-, high-mass) systems 
nor Roche Lobe Overflow (RLOF) binaries can explain the 
entire faint population. However, the quiescent transients 
may be responsible for most of the faint soft GC sources. 
Muno et al. (2004) suggested that the observed faint hard 
sources are most likely intermediate polars (IPs); a sub- 
class of magnetic cataclysmic variables (CVs). 

Among magnetic CVs, IPs are asynchronous rotators, 
with the white dwarf (WD) spin period being shorter than 
the orbital period. Currently there are ~ 31 IPs known 
(Gansicke et al. 2005). They consist of a magnetized WD 
and a low-mass companion. The companion transfers mass 
via RLOF to the WD, with typical mass transfer (MT) 
rates of ~ 10~ n M Q /yr. Matter spirals in toward the 
WD forming an accretion disc, which is truncated by the 
WD magnetosphere. Matter is then channeled into accre- 
tion columns over both magnetic poles of the WD (e.g., 
Belle et al. 2005). Typical IP magnetic fields are < 10 
MG (de Martino et al. 2004). It is from below the accre- 

1 Foreground sources excluded; M.Muno, private communication. 



tion shock, above the WD surface, where the hard X-rays 
originate (see Warner 1995; Patterson 1994 for review). 
IPs are known to exhibit both soft and hard X-ray emis- 
sion and are thought to be the most luminous subclass of 
CVs in X-rays, owing to their typically higher MT rates. 

We address the issue of the nature of the faint hard 
source population observed in the deep GC exposure of 
MM03, and test the validity of the Muno et al. (2004) hy- 
pothesis that these sources are IPs. We construct a simple 
phenomenological model of an IP and use population syn- 
thesis (§ 2) to calculate the number of IPs and their X-ray 
luminosities in the GC (§3). In §4 we comment on other 
population synthesis studies and discuss our results in con- 
text of available observations. 

2. MODEL DESCRIPTION 

We use the updated population synthesis code 
StarTrack (Belczynski, Kalogera & Bulik 2002; Belczyn- 
ski et al. 2005a). All stars are evolved with solar metallic- 
ity (Z=0.02) and with standard wind mass loss rates. We 
assume a continuous star formation rate in the GC over 
the last 10 Gyrs, and a binary fraction of 50%. Several 
physical processes important for binary evolution are ac- 
counted for: tidal interactions, detailed MT calculations, 
common envelope (CE) events, supernovae explosions, and 
various mechanisms for angular momentum losses such as 
magnetic braking (MB) and gravitational radiation (GR), 
among others. 

The goal of this study is to test the hypothesis that 
IPs can account for the observed population of the faint 
hard X-ray point sources in the GC. We use our standard 
model (Belczynski et al. 2005a) to check whether within 
the general framework of binary evolution, the number of 
predicted IPs coincides with that of the GC faint hard 
sources. It has been realized that CVs are the likely out- 
come of the CE phase (Paczynski, 1976). Despite years of 
investigation, no consistent physical model for this phase 
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of binary evolution exists. Therefore, we use two current 
alternative pictures of CE evolution to provide an estimate 
of associated model uncertainties. Our standard CE model 
calculations assume energy balance (Webbink 1984); the 
donor envelope is ejected from the system at the expense 
of the binary orbital energy. The binary separation fol- 
lowing the CE phase depends on parameters a cc and A, 
relating to the efficiency with which the orbital energy is 
used to expel the CE from the system, and donor internal 
structure, respectively. We adopt a cc x A = 1.0. In our 
alternative CE model, we have adopted the prescription of 
Nelemans & Tout (2005) which employs angular momen- 
tum balance, and assumes that the angular momentum is 
lost from the binary in a linear fashion as a function of 
mass loss. Following Nelemans & Tout (2005), we adopt 
a scaling factor 7 = 1.5 for this model. See Belczynski, 
Bulik & Ruiter (2005b) for the CE equations. 

X-ray Luminosity Calculation. We assume that the IP 
X-ray luminosity is a function of the accretion rate, accre- 
tor physical properties, and the efficiency with which the 
accretion luminosity is converted to hard X-ray luminosity 
in the Chandra band. 

For degenerate donors (M(i on ), we assume that MT is 
GR-driven and calculate it via 



M don = M don D- 
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where J or b is the orbital angular momentum, rfj, 
the angular momentum loss due to GR, / a is the 
tion of transferred mass accreted by the WD of mass 
M^cc (we assume Eddington limited accretion; see below), 
q = (M acc /M don ), and f3 mt = M acc M 2 on /(M don + M acc ) 2 . 
The radius mass exponent for the donor £ don is obtained 
from stellar models (see Belczynski et al. 2005a). 

For non-degenerate donors we estimate the MT rate as 
follows: 
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where £ C vi is the change of the donor radius due to its 
nuclear evolution, (j OD is the radius mass exponent for 
the donor Roche lobe, and T mD , r t ;d and the 
timescales associated with MB, tidal interactions and 
GR, respectively. In some cases MT proceeds on a 
thermal timcscale, and thus we use M t h = — M don T t h 
where the thermal timescale may be obtained from 
Tth = (30 x M d 2 on )/(i? don L don ) (e.g., Kalogera & Web- 
bink 1996). Then, we can calculate mass accretion rate 
as M acc = min(0.95 x M don , M odd ). We limit the accretion 
rate to the critical Eddington rate (M c dd), such that if 
MT is supcr-Eddington, the excess material leaves the bi- 
nary with the specific angular momentum of the accretor. 
We impose a mass-loss rate of 5% in the sub-Eddington 
regime since some IPs are observed to experience a mass 
loss rate of few percent (~ 10% in the extreme case of AE 



The hard Chandra band X-ray luminosity of an IP sys- 
tem is calculated from 
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where G is the gravitational constant, e is the conversion 
efficiency of gravitational binding energy to radiation (1 for 
WD surface accretion), and M acc and R acc are the accret- 
ing WD mass and radius, respectively. We use ?7boi = 0.09 
for the bolometric correction to 2-8 keV X-ray luminos- 
ity, although it is noted that this value is quite uncertain 
and may span a wide range (~ 0.01-0.2; Norton & Wat- 
son 1989b). The IP X-ray emission is likely anisotropic 
to some extent, since the accretion proceeds through the 
WD magnetic poles. However, it has been demonstrated 
that the emission region in IPs may be quite extensive 
(i.e., may encompass more than a quarter of the WD sur- 
face, Norton & Watson 1989a) and thus the anisotropy is 
not expected to be large. Here we assume isotropic emis- 
sion and so 77 gco = 1, although we note that a more so- 
phisticated model including polar cap accretion should be 
invoked once more observations are available. In general 
?7g o depends on the size of the emitting regions and their 
relative orientation to the observer. 

Initial Conditions. The initial mass of single stars and 
the primary (more massive) components in binaries are 
drawn within the range 0.8-150 M© from a broken three- 
component power law initial mass function (IMF), with a 
slope of-1.3/-2.2/-2.7in mass ranges 0.08-0.5/0.5-1/1-150 
M Q (Kroupa, Tout & Gilmore 1993). The secondary mass 
is obtained through a flat mass ratio distribution. Initial 
binary orbits are specified by semi-major axis (distribu- 
tion flat in the logarithm <~ 1/a) and eccentricity (thermal 
distribution ~ 2e). 

The stellar population of the GC is evolved through 
10 Gyr (age of the Galaxy) with a constant star for- 
mation rate and IP systems are extracted and luminosi- 
ties compared against the MM03 survey sources. We 
call any binary system experiencing RLOF in which the 
accretor is a WD and the donor is any type of star a 
CV (i.e. including AM CVn systems, see Warner 1995). 
We consider various types of WDs: helium (HcWD), 
carbon-oxygen (COWD), oxygen- neon (ONcWD), hydro- 
gen (HWD; formed by stripping the envelope of a low-mass 
main sequence (MS) star in RLOF), and hybrid (HybWD; 
containing a He envelope and a He-C-0 mantle, formed by 
stripping of a low-mass He star in RLOF). We note that 
the evolution of IPs and CVs are identical in our model, 
i.e., we do not incorporate the magnetic field of the WD 
into the evolution. Rather, from the CV population we as- 
sign a fraction of CV systems that are IPs. In our model 
we adopt an IP fraction (IP frac ) of 5% following MM03. 
We have calibrated our results to pertain to the surveyed 
GC region by scaling by stellar mass. Following Muno et 
al. (2004; see their §3), we have assumed that the 17'xl7' 
Chandra field of view of the GC corresponds to a cylinder 
440 pc deep with a radius of 20 pc, encompassing 1.3 x 10 8 
M Q in stars 2 . 

3. RESULTS 



Aqr; Wynn, King & Home 1997). 

2 There may be 50% uncertainty associated with this estimate (Launhardt et al. 2002), which will propagate linearly in our results. 
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It is found that the synthetic GC population of IP sys- 
tems span a range of X-ray luminosities ~ 3x 10 29 — 5x 10 33 
erg s _1 . In the following we only discuss systems above 
the MM03 survey detection limit, i.e., with X-ray lumi- 
nosities > 10 31 erg s _1 unless otherwise noted. The main 
results of our calculations for the standard and alternative 
CE prescriptions are presented in Tableland Figure [I] 

Standard CE model. The number of IPs depends 
strongly on the adopted IP frac . We find - 800-8000 IPs 
in the GC for IP nac of 1-10%, respectively. To match the 
observed number of faint hard X-ray sources in the GC 
(1427) with IPs we would require an IPf rac of ~ 2%. The 
distribution of orbital periods of our IPs peaks between 
~ 1.5-3.5 hr and extends toward longer periods (> 5 hr), 
which is in agreement with the observed orbital periods of 
magnetic CVs (see i.e., Tovmassian et al. 2004). We find 
that the most frequent IP type (89%) is a magnetic WD 
accreting from a MS star. The majority of donors (87%) 
are < 0.7 M s , i.e. K-M dwarfs. Accretors in the WD-MS 
subclass are COWDs (66%) and HeWDs (23%). The sec- 
ond most frequent (10.5%) type of IP is a double degener- 
ate system (i.e., close WD- WD binary with stable RLOF). 
The most common configurations found within this sub- 
class are: COWD-HeWD (6.6%), COWD-HybWD (1.7%) 
and COWD-HWD (1.7%). Finally, we find few (0.5%) 
IPs with giant- like donors. In this subclass nearly all sys- 
tems consist of a low-mass (~ 0.3 — 0.4 M Q ) sub-giant or 
a red giant transferring mass to a COWD. The relative 
occurrence frequencies of different IP types reflect various 
evolutionary timescales for donors of varying types and 
masses. The WD-MS IPs are abundant since low-mass 
MS stars have very long evolutionary timescales (~ 10 11 - 
10 12 yr) and RLOF at the IP phase is also driven on long 
timescales; by GR (rg r ~ 10 10 yr) with the addition of 
MB (r m b ;> 10 9 — 10 yr) for some systems. The evo- 
lution of double degenerates is driven by GR on shorter 
timescales on the order r gr ~ 10 9 yr (since these binaries 
are tighter); moreover it is difficult to form these systems 
(survival of two CE phases rather than one). For WD- 
giant donor systems, the evolutionary and MB timescales 
are shorter than in the above cases, making them the least 
represented subclass of IPs. 

In Figure ^ we show the overall luminosity distribution 
and the corresponding X-ray luminosity function (XLF) 
for the model IPs. The IPs above the MM03 detection 
limit are marked for easy comparison. It is noted that 
only ~ 55% of IPs are bright enough to make the X-ray 
luminosity threshold of the survey. The brightest IPs in 
our simulations are those with giant donors, with average 
X-ray luminosities L x ~ 5 x 10 32 erg s _1 . The power-law 
slope of the cumulative distribution (N(> L x ) ~ L~@) of 
the XLF for the synthetic IPs (above the survey detection 
limit) is p ~ 0.8. This value is shallower than the one 
derived by MM03 3 (1.34 for L x > 4 x 10 31 erg s" 1 ) and 
Muno et al. (2006). It is noted that the full XLF slope 
modeling for IPs would require i) taking into account ab- 
sorption toward the GC sources (as well as circumbinary 
absorption) , and ii) a detailed model of anisotropic X-ray 
emission from a magnetized WD. 

Alternative CE model. It is found that ~ 170-1700 IPs 

3 Although formula (5) in MM03 does not represent a cumulative 
communication. 



may be present in the GC for IPf rac of 1-10%, and one 
would require an IPf rac of ~ 8% to match the number of 
GC faint hard sources. Orbital periods are found in the 
same range as for the standard CE model. Once again, the 
most frequent IP type is a magnetic WD accreting from a 
low-mass MS star (90.8%) with the majority of them be- 
ing COWD-MS star binaries (75%). Double degenerates 
are the second most populated subclass (8.6%), first with 
COWD-HWD (4%), followed by COWD-HeWD and then 
COWD-HybWD (2.1% and 1.9%, respectively) with other 
types constituting the rest. Again only a small fraction 
(0.6%) of the IPs involve a WD-giant star binary. The 
slope of the XLF for IPs above the survey detection limit 
is found again to be P ~ 0.8 (see Fig 1). 

Additionally, we have calculated a model with decreased 
CE efficiency (a ce x A = 0.1). The number of IPs decreases 
by ~ 50%, due to the more frequent binary mergers during 
CE phases. Also, we have tested the assumption on the 
initial mass ratio of binary components with a model in 
which both binary component masses are chosen indepen- 
dently from the IMF. The overall binary evolution is quite 
different, since on average q is much smaller (~ 0.1, as op- 
posed to the flat g-distribution in other models). There is 
an overall increase (~ 20%) in the number of IPs, though 
the population content is very similar with 89% of IPs with 
> L x 10 31 erg s _1 being MS- WD binaries. 

Typical evolution. Two MS stars (2 and 0.6 Mq) start 
out on a highly eccentric orbit (e=0.9) with an orbital pe- 
riod of 1460 days when the Galaxy is 2.1 Gyr old. At 3.6 
Gyr the system has circularized (period of 109 days), a CE 
phase is initiated by the primary (now an asymptotic gi- 
ant), and upon ejection of the envelope the orbital period 
shrinks by ~ two orders of magnitude (now 16 hours). The 
primary shortly thereafter becomes a COWD. At 8.4 Gyr 
the MS secondary initiates RLOF and the system becomes 
an IP. At 10 Gyr we find an IP system with a period of 2.9 
hours. We note that for this same system, if the alterna- 
tive CE prescription model is used an IP is never formed. 
The reason for this is that upon ejection of the CE, the 
binary does not lose enough orbital angular momentum in 
order to end up in a tight enough orbit such that RLOF, 
and an IP phase, can ensue. Hence a smaller number of 
IPs are found in our alternative CE model. 

4. DISCUSSION 

We have explored the possibility that the faint hard 
sources in the GC are IPs. It is found that for both current 
CE models, the IP population is ample enough to explain 
the GC faint hard sources. The required IP fractions are 
then ~ 2% and ~ 8% for the standard and alternative CE 
models, respectively. The actual fraction is uncertain and 
currently estimated to be ~ 5% by MM03 based on the 
Kube et al. (2003) catalog of CVs. Once more stringent 
observational constraints on the IP fraction are obtained, 
we will be able to test the validity of different CE models. 

At this point a full picture of the GC X-ray point sources 
begins to emerge. Most of these sources are faint and are 
spectrally hard (1427) and can be explained by a popu- 
lation of IPs as suggested by Muno et al.(2004) and con- 
firmed in this study. Additionally, there may be a small 

ibution, the listed value of the slope is cumulative; M.Muno private 
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contribution (~ few percent) of wind-fed sources with NS 
and BH accretors to the faint hard population, as proposed 
by Pfahl et al. (2002) and later revised by BT04 and Liu 
& Li (2005). The faint soft sources (652) are likely RLOF 
transients with NS and BH accretors in quiescence as pro- 
posed by BT04. The bright GC sources 20; Wang et 
al. 2002 and MM03) can be explained by a population of 
NS/BH persistent sources and transients in outburst (e.g., 
BT04). The stellar density of the GC region is higher 
than that of the solar neighbourhood, so it's possible that 
dynamical interactions may alter the number of CVs in 
the GC. Tidal captures or exchange interactions may in- 
crease the numbers, but also, the primordial CV popula- 
tion may be depleted by tidal disruptions in dense envi- 
ronments (i.e., see Heinke et al. 2003). A full evolutionary 
model of the GC incorporating binary stellar evolution and 
dynamical interactions should be invoked to gauge these 
processes. Here, we have just demonstrated that field-like 
(no stellar interactions) evolution provides enough IPs to 
explain the GC faint hard X-ray population. 

Liu & Li (2005) conducted a different population synthe- 
sis study of the GC sources, and it was found that most 
of the Chandra sources are NS Low Mass X-ray Binary 
transients, with IPs playing a minor role in the Wang et 
al. and the MM03 surveys. The IP MT rates obtained by 
Liu & Li are much lower (~ 2 orders of magnitude) than 
our predictions and those of observational estimates. This 
leads to an underestimate of IP X-ray luminosity in the 
Liu & Li (2005) simulations, and hence a severe underesti- 
mate of IP numbers in the GC. The Liu & Li (2005) model 



is based on the population synthesis code by Hurley, Tout 
& Pols (2002). MT is calculated from a formula that does 
not directly incorporate the star's response to mass loss, 
nor the various mechanisms for angular momentum loss. 
It was originally obtained and calibrated for binary sys- 
tems with evolved donors (Algols; Tout & Eggleton 1988) 
and does not appear to work for low mass IP donors. 

We have found that most of the Galactic Centre IPs 
are either magnetic WDs feeding from low-mass MS late- 
type companions, or double degenerate systems (e.g., AM 
CVn). These systems are very unlikely to be detected at 
wavelengths other than X-rays due to the high extinction 
toward the GC, and the low luminosity of the IP systems 
(K-M dwarfs). Recently, Laycock et al. (2005) carried 
out a search for infrared counterparts of the GC X-ray 
sources. It was found that high mass X-ray binaries (with 
donors brighter than B2V) are unlikely candidates for the 
majority of the GC faint sources, in agreement with our 
findings and those of BT04. Bandyopadhyay et al. (2005) 
are conducting a deeper near-infrared survey which will 
detect all giant type donors and MS donors with spectral 
types earlier than G. However, if the typical donors in IPs 
are indeed K-M MS stars, as found in our study, and also 
observed for most Galactic IPs (e.g., M3 V for EX Hya, 
Dhillon et al. 1997; K5 V for AE Aqr, Tanzi et al. 1981), 
these systems will go undetected in this survey. 

We acknowledge the support of KBN grant 
1P03D02228, and thank M.Muno, J.Grindlay, K.Belle, 
X.-W. Liu & X.-D. Li for helpful discussions. We also 
thank the anonymous referee for insightful comments. 
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Table 1 

Galactic Centre Population of IPs 



Model 


WD-MS 


WD-giant 


WD- WD 


Total 


Standard CE 










IPfrac 1% 


694 


5 


82 


781 


IPfrac 5% 


3474 


25 


410 


3909 


IP frac 10% 


6948 


50 


820 


7818 


Alternative CE 










IPfrac 1% 


157 


1 


14 


172 


IPfrac 5% 


784 


5 


74 


863 


IPfrac 10% 


1568 


9 


149 


1726 




Fig. 1. — Left: X-ray luminosity distribution and corresponding XLF for GC IPs for the standard CE model and an IP 
fraction of 5%. The entire population is shown with a dotted line while IPs brighter than the MM03 X-ray luminosity 
threshold are shown with a solid line. The XLF slope of the cumulative distribution (N(> L x ) ~ L~P) is marked. Right: 
same as the left panels but for the alternative CE model. 



